Quantum plasmonics with a metal nanoparticle array 
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We investigate an array of metal nanoparticles as a channel for nanophotonic quantum communication and 
the generation of quantum plasmonic interference. We consider the transfer of quantum states, including single- 
qubits as plasmonic wavepackets, and highlight the necessity of a quantum mechanical description by comparing 
the predictions of quantum theory with those of classical electromagnetic theory. The effects of loss in the metal 
are included, thus putting our investigation into a practical setting and enabling the quantification of the perfor- 
mance of realistic nanoparticle arrays as plasmonic quantum channels. We explore the interference of single 
plasmons, finding nonlinear absorption effects associated with the quantum properties of the plasmon excita- 
tions. This work highlights the benefits and drawbacks of using nanophotonic periodic systems for quantum 
plasmonic applications, such as quantum communication, and the generation of quantum interference. 
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I. INTRODUCTION 



The field of quantum plasmonics is currently experiencing 
intense interest from the plasmonics and quantum optics com- 
munities 1 1-24]. Integrated quantum systems featuring sur- 
face plasmons are showing remarkable potential for their use 
in quantum control applications, such as quantum informa- 
tion processing [12-16]. Here, novel capabilities in the way 
the electroma gne tic field can be localized [|25l [26ll and ma- 
nipulated [|27l - [3lll offer the prospect of miniaturization, scal- 
ability and strong coherent coupling to single emitter sys- 
tems that conventional photonics cannot achieve lfT6l - l24l]. Re- 
cent studies have focused on entanglement preservation | US], 
quadrature-squeezed surface plasmon propagation ^ and the 
use of surface plasmons as mediators of entanglement be- 
tween two qubits [4,-^]. With the advancement of nanofab- 
rication techniques, ordered arrays of closely spaced noble 
metal nanoparticles have been proposed as a means of guid- 
ing electromagnetic energy, via localised surface plasmons 
(LSPs), on scales far below the diffraction limit ll32l[33ll . Here, 
energy transport relies on near-field coupling between surface 
plasmons of neighbouring particles [34], with the suppression 
of radiative scattering into the far-field li35i437[l . Recently 
it was shown that an appropriate arrangement of nanoparti- 
cles can form passive linear nanoscale optical devices such as 
beam splitters, phase shifters and crossover splitters if sslfioll . 
While much progress has been made in the area of device de- 
sign, so far there has been no analysis of the effects of loss in 
these nanoparticle systems in the quantum regime. It is vital 
to understand the impact of these effects on the performance 
of such devices so that plasmonic systems may be developed 
as an efficient platform for nanophotonic quantum control ap- 
plications. 

In this work we carry out such an analysis and investigate 
quantum state transfer and interference of surface plasmons 
on a metal nanoparticle aiTay. The transfer of quantum states, 
including those encoded into single-qubit plasmon wavepack- 
ets, is studied. The effects of loss in the metal due to elec- 
tronic relaxation are also included in our model, putting the 
investigation into a more practical setting. We find that quan- 



tum state transfer can be achieved for small length arrays even 
under nonideal conditions and therefore these arrays may act 
as channels for short distance on-chip nanophotonic quantum 
communication. We also study the interference of single plas- 
mons in the nanoparticle array and find nonlinear absorption 
effects associated with the quantum properties of the plasmon 
excitations. Our study highlights the benefits and drawbacks 
associated with building nanophotonic systems that use sur- 
face plasmons in the quantum regime. The results of this work 
may help in the future study and design of more complex plas- 
monic structures involving emitter systems for quantum con- 
trol applications, and the probing of novel nanoscale optical 
phenomena. 

We start our investigation in the next section by introduc- 
ing the nanoparticle aiTay model and quantized mathematical 
description, along with some basic properties of the system 
dynamics. Then in section III we study the performance of 
quantum state transfer under ideal conditions and highlight the 
necessity of a quantum mechanical description by comparing 
the predictions of quantum theory with those of classical elec- 
tromagnetic theory. In Section IV we consider the effects of 
damping due to losses associated with the electronic response 
within the metallic nanoparticles and study the interference of 
single plasmons in the nanoparticle array. Finally, we summa- 
rize our findings in Section V. 



II. PHYSICAL SYSTEM AND MODEL 

A. The Hamiltonian 

We consider the system depicted in Fig.[T](a), which is pre- 
sented in a top-down view. Here, a tapered metal nanowire 
waveguide on the left hand side focuses light at the end of 
its tip in the form of a confined surface plasmon field. This 
field then couples to the adjacent spherical metal nanoparticle 
and excites a localised surface plasmon (LSP). The LSP ex- 
citation propagates across the linear array of metal nanoparti- 
cles by near-field coupling and exits via another tapered metal 
nanowire waveguide on the right hand side. All metal re- 
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FIG. 1: (a): A tapered metal nanowire waveguide on the left hand side focuses light to the end of its tip and excites a localised surface 
plasmon (LSP) on the adjacent nanoparticle. The excitation then propagates across the array of nanoparticles and exits via another tapered 
metal nanowire waveguide on the right hand side. All metal regions have permittivity e,„ and dielectric regions have permittivity 6d, as defined 
in the text. The dimensions are chosen as an example and the theory developed is more general, with a range of parameters investigated in 
this work, (b): Weak coupling approximation for the nanoparticle array. The theoretical model developed in this paper is valid in the regime 
1^1 <IC Wo and in particular we choose max|g| = O.loio, where coo is the natural frequency of the nanoparticle field oscillations and g is the 
nearest neighbour coupling parameter. The lower blue curve is for transverse polarization (T) and upper red curve is for longitudinal (L) 
polarization. 



gions have a frequency dependent permittivity e,„(w) and di- 
electric regions have static real and positive permittivity ej. 
In Fig. [T] (a), we give a specific example of the system be- 
ing studied by choosing the radius of the nanoparticles as 
R = 25nm and the distance between nanoparticles in the ar- 
ray as c/ = 75nm (however, the general model we will intro- 
duce allows arbitrary values to be chosen for these parameters 
and for all other physical parameters). We consider the metal 
nanoparticles in the array support electron charge density os- 
cillations in the longitudinal (L) and transverse (T) directions, 
as shown in the inset of Fig. [T](a) and neglect multipolar inter- 
actions ll4lll . In the system depicted in the main part of the fig- 
ure, due to the direction in which the electron charge density 
oscillates in the nanowires, which act as the surface plasmon 
source and drain on the left and right hand sides, the nanotips 
at the ends are oriented to excite/collect charge density oscil- 
lations in the L direction. For excitation and collection in the 
T direction, both nanowires should be rotated by 90 degrees 
either clockwise or anticlockwise. Further details regarding 
the nanowire orientations are discussed later. 

We now introduce the Hamiltonian for the system, justify- 
ing the physical origin of each of the terms appearing. The 
total Hamiltonian describing the system in Fig.[T](a) is given 
by 

H - H„p + Hs + Hd + H„p^s + H„pj. (1) 

Here, the first term describes the linear nanoparticle array 
system consisting of n nanoparticles and is given by 

n 

H„p = ^ htojajai + ^ hgij (fl|a j + aja,) , (2) 
i=i [ij] 

where is the natural frequency of the field oscillation at 
the i-th nanoparticle, gij is the coupling strength between the 



fields of the /-th and j-th nanoparticles, [/, j] denotes a sum- 
mation over nearest neighbours j for a given nanoparticle 
and the operators a! (a,) represent the creation (annihilation) 
operators associated with a field excitation at nanoparticle site 
/ which obey bosonic commutation relations [a,, aj] = 6ij. 
Here, a macroscopic quantization of the fields is used, where 
the field modes are defined as localized solutions to Maxwell's 
equations satisfying the boundary conditions of the metal- 
dielectric interface (4^. In this case, the electron response 
is contained within the dielectric function of the metal SEH]. 
We consider either L or T polarization along the array, sup- 
pressing the polarization index. In addition, while the model 
we investigate here is for a linear array of nanoparticles, the 
theory introduced can be applied to more complex arrange- 
ments of nanoparticles [38-40]. 

The first term in Eq. (|2|l represents the free Hamiltonian of 
the fields at the nanoparticles, where satisfies the Frohlich 
criterion, Re[e,„(dLt,)] = ^2ed [34, 37]. This criterion con- 
siders the nanoparticles to be small enough compared to the 
operating wavelength such that only dipole-active excitations 
are important ll27ll . Taking all nanoparticles to have the same 
permittivity e,„, the local frequencies can be set to be equal, 
w, - cl>[), V/. Due to the spherical symmetry of the nanopar- 
ticles, these local frequencies are independent of the polar- 
ization. The second term in Eq. (|2]l represents a nearest- 
neighbour coupling between the near-field at each nanoparti- 
cle. In order to justify the physical mechanism of this second 
term, we briefly provide the correspondence of the quantum 
description of the nanoparticle array to the classical descrip- 
tion Q. 

Consider a quantum state \ifr) - Yli where la,) = 
e^^l"'l"e"'" |0) is a coherent state and a,- is the mean field 
amplitude at the i-th nanoparticle. Here, the electric field vari- 
ation of a coherent state \a) approaches that of the classical 
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wave picture in the Umit of large amplitude a 1*47] . Taking 
H„p and and substituting them into the Schrodinger equa- 
tion, ihj^ \if)) - H„p It//), one finds the differential equation for 
the mean field amplitudes as llssll 



dai 
dt 



[ij] 



(3) 



By choosing all the couplings to be equal gi j - g - 
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where y - Jt - ^ and y - ji - —2 are the relative 
couplings and phases for polarization T and L respectively 
(at a fixed distance d, array orientation and nanoparticle size 
R lH^l), the differential equation in Eq. (O is exactly the same 
as the classical differential equation for the amplitude of the 
dipole moment p, (associated with the electric field at site /) 
for an array of interacting Hertzian dipoles under the condi- 
tion 0)1 ^ Wo for the interaction frequency w/ |34, 38t]. This 
is a weak coupling approximation including only the nearest- 
neighbour interactions. In the classical Hertzian model, the 
dominant interaction in the system is considered to be between 
the nanoparticle dipoles via the Forster field, which has a l/t/-* 
dependence for d <^ A, where A - Aq/ ^fea and Aq is the free- 
space wavelength corresponding to the natural frequency o)q 
of the nanoparticle dipole field. An — IticIojq (c is the ve- 
locity of light in a vacuum) fs^BH. This regime {d <C /I) 
is known as the near-field approximation. Furthermore, the 
dipoles are considered point-like for R < dl3 |41], known as 
the point-dipole approximation. Thus, under the weak cou- 
pling, near-field, and point-dipole approximations, the quan- 
tum model with g - o)jy/2o)Q recovers the classical dynamics 
in the correct limit using coherent states. Here, the interac- 
tion frequency is given by ojj — [e^p^iR^ /3rn* eoedd^]^^^ i34ll . 
where e is the electronic charge, pei is the free electron density 
of the metal, m* is the optical effective electron mass and eo is 
the free-space permittivity. 

In Fig. [1] (b) we show an example of the dependence of 
the magnitude of the coupling g = a;^y/2wo (in units of 
o)o) as the ratio of R/d increases. Here we have taken the 
permittivity of the metal e,„(w) as silver and used €,nio)) - 
foo - Wp/(w^ + iwT) + ill, where = 5, F = 6.25 x lO'^ 
rad/s and o)p - 1 .402 x 10"' rad/s, which are chosen to obtain 
a best fit to experimental data at frequencies corresponding 
to freespace wavelengths Aq > 350nm ll46ll . i.e. the optical 
range and above. This leads to wq = 5 x lO'^ rad/s, the lo- 
cal frequency of the nanoparticles. In addition, we have used 
Pel = 5.85 X lO^'^m-^ m* = 8.7 x lO^^'kg and ea = 1 IH. 
The weak coupling approximation is equivalent to |g| ^ wo 
and we impose this by setting max|g| - O.Iwq. Note from 
Fig. [T] (b) that the condition max|g| = O.lwo satisfies the 
point-dipole approximation immediately, as well as the weak 
coupling for both polarizations. For the near-field approxima- 
tion to also be satisfied we require d ^ Ittc/ojq ~ 377nm. 
The example in Fig.[T](a) with d - 75nm and R - 25nm with 
silver satisfies all three of the required approximations. 

An additional requirement for the system is that quantum 
effects other than those due to the quantized surface plasmon 
field, such as electron tunneling between nanoparticles and 
the quantum size effect of each nanoparticle [47], are neg- 
ligible. This puts a lower limit on the distance d between 
nanoparticles at ^ Inm ll48ll . and nanoparticle radii of the or- 



der of Inm (37^, respectively. However, in order to confidently 
use the macroscopic approach for the quantization of the sur- 
face plasmon field due to the electron response, we assume 
nanoparticle radii R > lOnm and therefore d > 30nm to sat- 
isfy the point-dipole approximation. As far as we are aware it 
is still an open question as to what dimension the macroscopic 
approach to surface plasmon quantization breaks down. In ad- 
dition, for the moment, we also neglect internal electronic re- 
laxation at the nanoparticles and relaxation of the dipoles into 
the far-field. Damping will be introduced after the ideal case 
has been developed in the next section. 

Continuing with our description of the physical system, the 
second and third terms of Eq. ([T]i represent the free Hamil- 
tonian of surface plasmon fields in the source and drain 
nanowires on the left and right hand sides of Fig. [T] (a) re- 
spectively and are given by 

Xoo 
do)fia>s\Lo)s(o)), 
-00 

Xoo 
do)tia>d^ (o))d(a>) . 
-00 

The operators of the nanowires correspond to continuum 
modes of surface plasmons, which obey the bosonic commu- 
tation relations [s(o)), s^{o)')] = 6{cL>—a)') and [d(a)),d\cL>')] = 
6{cl> — to'). Again, a macroscopic quantization is carried out 
for the field Q [23l 14311 and we have extended the integration 
of 0) to cover the range —00 to 00 ll44ll . 

The surface plasmon excitation in each nanowire is taken 
to correspond to the fundamental transverse magnetic mode 
with winding number m - |23]. It can be generated by 
various methods. For instance, it could be generated via cou- 
pling of a photon from the far-field by focusing the quantized 
light field onto a grating structure at a thicker part of the ta- 
pered nanowire |4y] . Another method could be to use end-fire 
coupling of photons in conventional silica waveguides to the 
metal nanowires, again at a much thicker part of the tapered 
wire 1 50]. One could also generate the plasmon excitations di- 
rectly on the wires very close to the tip region by driving emit- 
ter systems, such as quantum dots [17] or NV-centers II19(] . 
with coherent light to further reduce losses during propaga- 
tion of the input ||23|] . Here, the combination of metal and 
emitter system may provide additional flexibility in optimiz- 
ing the field profile that couples to the nanoparticle - similar 
to a nanoantenna system ifsTl] - rather than a direct coupling of 
the emitter system on its own. 

The fourth and fifth terms of Eq. ([T]i represent the coupling 
of the surface plasmon field of the source nanowire to the LSP 
field of nanoparticle 1 and the surface plasmon field of the 
drain nanowire to the LSP field of nanoparticle n respectively. 
Using a weak-field linearized model 1521 153[] the terms are 
given by 



Hnp,s = ifi 



f 



da)gin((^)[a\s' (lo) — s{o))a'.], 



Xoo 
da)gouiio))[a„d^ (o)) - d(a))al]. 
-00 

Here the coupling parameters gin/out(w) depend on the strength 
of the near-field coupling between the nanowires and nanopar- 
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tides. Focusing on the case of the source nanowire-to- 
nanoparticle coupling and considering a propagating surface 
plasmon in the nanowire entering the region at the nanotip 
from the left hand side, we have that for an appropriate 
paraboloidal profile of the nanowire, the mode function of the 
excitation near the tip strongly couples to a dipole orientated 
in the direction of the propagation and placed in close prox- 
imity |23, 54]. Thus, for the orientation of the source nanotip 
shown in Fig.[T](a), the nanowire field couples predominantly 
to the L polarized oscillation in the nanoparticle. For coupling 
to the transverse polarization, we rotate the nanotip clockwise 
by 90 degrees. The reciprocal case holds at the drain nan- 
otip and in the orientation shown in Fig. [T](a), the drain pre- 
dominantly couples to L polarized field oscillations in the «-th 
nanoparticle. 

Regardless of the excitation method of the plasmons in the 
nanowire, for a wire that has a slowly varying radius R{x) with 
distance x from the tip along the wire, we have a locally vary- 
ing dispersion relation given by fssll 

em Ii(kQK„R(x)) ^ Q Ki(koKdR(x)) _ ^ 

Km Mh)KmR{x)) Kd KoikoKdRix)) 

where Ip and Kp are the modified Bessel functions, - 

-ijn^ — e,„, Kd - s/n^ — Cd, = w/c is the freespace 
wavenumber at a given frequency co and n - n(x, at) is the lo- 
cal effective refractive index at position x along the nanowire. 
The radius of the nanotip at the end of the source wire de- 
fines the effective radius of the wire in the region just before 
the tip. We can use this to determine the approximate dis- 
persion relation of the surface plasmons entering the nanotip 
region where they couple to the LSP of the first nanoparticle. 
Therefore, Eq. (|4]l can be solved for a given set of physical 
parameters in order to obtain the local effective refractive in- 
dex n(x, (jj), leading to the dispersion relation k - n{x, w)w/c 
for the surface plasmons close to the tip (x ~ 0). In Fig. |2] 
we show the dispersion relation for a free-space photon k - 
k() = a>/c, a nanowire surface plasmon (for R(Q) = 25nm) 
using Eq. (|4]l, a standard metal-air interface surface plasmon 
k = (oj/c) \/e„{(L))/{l + e,„(w)) |8] and the nanoparticle natu- 
ral oscillation frequency ojo- In all cases the example metal 
is taken to be silver with the dielectric function defined previ- 
ously. Here we have chosen to represent the wavenumber k in 
units of an array spacing d = 75nm. Note that only the sur- 
face plasmon field from the tip region of a nanowire has the 
potential to achieve both the correct energy conservation (co 
matching) and dipole-coupling ll23ll54|] for eflicient near-field 
coupling to the first (or last) nanoparticle of the array. 

Thus, by setting the coupling gmioS) in H„p^s according to 
the physical geometries and tip orientation being considered, 
one can model coupling of the surface plasmon in the source 
nanowire to the first nanoparticle and its reflection back along 
the nanowire. Similarly, by setting the coupling gouti^^) in 
Hnp^d, one can model coupling of the last nanoparticle's near- 
field to the drain nanowire and its reflection back along the 
nanoparticle array. As the field profiles at the tips are sim- 
ilar in form to those of the nanoparticles, the same physi- 
cal approximations as the inter-particle coupling strengths gij 
should be satisfied by the couplings ^in(w) and gom((^) in or- 
der for the model to be a consistent description. The gin/out(w) 
couphngs can then be modified to model non-ideal mode func- 
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FIG. 2: Dispersion relation for various forms of surface plasmon 
excitation. Here silver has been chosen as the metal and air as the 
dielectric background media. The straight dashed line P shows the 
photon dispersion relation in free-space, the curve SPP-f is the dis- 
persion relation of a surface plasmon field at a standard metal-air 
interface and the curve SPP-w is the dispersion relation of a surface 
plasmon field in the tapered metal waveguide in the region of its tip 
at a radius of R = 25nm. The horizontal dotted line np is the nat- 
ural frequency of a single nanoparticle a>o. For all curves only the 
real part of the wavenumber k is plotted. The imaginary part being 
several orders of magnitude smaller. 

tion profiles due to the tip shape and other geometrical factors. 



B. Transmission and dispersion 

We now use the Hamiltonian in Eq. ([TJ to model the trans- 
mission of a quantum state injected into the array by the 
source nanowire and then its propagation along the array until 
it is subsequently extracted out by the drain nanowire. In or- 
der to do this we use an effective scattering matrix approach 
that will link the input field operators of the source nanowire 
to the output field operators of the drain nanowire, provid- 
ing a method to map arbitrary input quantum states to output 
quantum states. This scattering matrix is obtained by apply- 
ing input-output formalism l52l Issll to the nanoparticle array, 
as summarized in Appendix A. The benefit of this approach is 
that we may treat the nanoparticle array as a waveguide with 
an effective medium, which makes the description of the sys- 
tem in the context of the transfer of quantum states more intu- 
itive. However, it is important to note that one can also use this 
approach to investigate the internal quantum dynamics of the 
nanoparticle array and even interactions with other resonant 
systems. For instance, emitter systems such as NV centres, 
placed in close proximity |52, 53]. 

Using the input-output formalism in Appendix A, we 
have the relation between input field operators Si„{cij) and 
iiin(a;), and output field operators iout(w) and douiicS) for the 
nanowires as follows, 

iin(w) = 7^*(w)^out(w) + 7;*(w)Jout(w), (5) 

dinio)) = 7^*(dL»)Sout('y) + 7?.^(w)Jout(w), (6) 

where the transmission Ts,d and reflection TZs^d coefficients are 
functions of the system parameters gi^, gout, gij and w,-, and 
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the relation |7?,5_^(a;)p + |Ts,rf(w)p - 1. Taking the Hermi- 
tian conjugate of Eq. (|5]l we have ij^Ciw) - Ti-s(<^)^l,M^'^^ + 
7^(w)(i^m(w). This allows us to describe the nanoparticle ar- 
ray as an effective waveguide, with transmission %(a)) - 
|7^(w)|e'^*''^'^', where x - (n + l)d is the total effective dis- 
tance (from the centre of the source tip to the centre of the 
drain tip, as shown in Fig. [T] (a)). The factor ztn takes into 
account the phase difference of — 1 in the definition between 
the input and output field operators and the wavenumber k de- 
pends on the system parameters ^in, gout, gij and w,-. We now 
drop the index s in the transmission coefficient for ease of no- 
tation and consider only transmission in the forward direction. 
Thus, with the use of T{cl>), we obtain 



arg [r(w)] T ^ + 2m7r n _li _lo 

k = , m = 0, ±1, ±2, 

(n + l)d 



(7) 



where the additional factor of Imn is included to reflect the 
cyclical degeneracy of the wavenumber. 

In Fig.[3](a)-(e) we plot the ampUtude squared of the trans- 
mission, |T(w)p, as the frequency is varied for an array of 
n = 1, 2, 3, 5 and 7 nanoparticles respectively. While 
an analytical form for T(tij) can be found, due to the gen- 
eral complexity of all the system parameters, here we show 
only explicit examples where we have taken all local frequen- 
cies to be equal iw, = (^o,yi and the couplings to be equal 
gij - gnp - — O.lwcV/ and its nearest neighbors j (with the 
minus sign for the longitudinal polarization, as yi = —2). In 
Appendix B we provide the analytical form for the T(w)'s. 
Physically, this chosen coupling regime corresponds to an ar- 
ray with d/R ~ 3, for example R = 25nm and d = 75nm, 
if we take the metal to be silver as before. The source and 
drain couplings are set as gin/out = 0.0 Iwo, achieved by vary- 
ing the distance between the nanowire tips and their respec- 
tive nearest nanoparticle. For a given number of nanoparticles 
n, the transmission spectral profiles in Fig. |3] (a)-(e) have n 
resonances at frequencies Luir^ = a>o + 2g„pC0s(k-'d), where 
y - jnlin +\)dfoTi i - \,...,n. In Fig.[3](f)-(j) we plot the 
effective wavenumber k from Eq. (|7]i for « = 1, 2, 3, 5 and 7 
respectively. Points corresponding to the y transmission res- 
onance peaks from Fig. [3](a)-(e) are marked as circles. Also 
included in these figures is the dispersion relation for the in- 
finite array case (dashed line), where the k^ take on continu- 
ous values [38i1, with a positive group velocity over the entire 
range due to taking the minus sign (phase) for the longitudi- 
nal polarization coupling g in this example. From Fig.[3](f)-(j) 
one can clearly see that as n is increased, the band structure 
of the infinite array case is gradually recovered, where each 
{(Orp k^) point corresponds to the dominant excitation of a sta- 
tionary eigenstate of the system Hamiltonian H„p, analogous 
to the case of coupled cavities lf56ll . for instance in the case of 
photonic crystals iIstIi . 

Note that while the above examples provide a basic in- 
sight into the system dynamics, the formalism introduced 
here can be used to describe more complex and general plas- 
monic nanoparticle systems with arbitrary couplings and nat- 
ural local frequencies. We now proceed to focus on odd num- 
bered nanoparticle systems with n > 1 in order to understand 
the transmission properties of larger arrays in more general 
regimes. A similar study could be made for even numbered 
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FIG. 3: Transmission spectral profiles and dispersion relations for 
undamped arrays of nanoparticles. Panels (a)-(e) correspond to the 
amplitude squared of the transmission, |T(<y)p, as the frequency is 
varied for an array of n = 1, 2, 3, 5 and 7 nanoparticles respectively. 
Panels (f)-(j) correspond to plots of the effective wavenumber k for 
;? = 1, 2, 3, 5 and 7 respectively. Also shown are points correspond- 
ing to the transmission resonance peaks from panels (a)-(e) as well 
as the dispersion relation for the infinite array case (dotted line). 



systems, however, we choose odd numbered as there is always 
a resonance at the natural frequency wq- This will become im- 
portant later in our study of quantum state transfer. 

In Fig.|4](a), (e) and (i) we plot a cross-section of the ampli- 
tude squared of the transmission, |T(w)p, as the frequency to 
and coupling gin (- gout) is varied for an array of « = 3, 5 and 
7 nanoparticles respectively. Here, Aoj - o) — uq and we have 
chosen to plot all parameters in units of the nanoparticle cou- 
pling gij - gnp, V/ and its nearest neighbors j. The plots are 
therefore independent of gnp, as long as gnp ^ cjq is satisfied. 
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FIG. 4: Cross-sections of the amplitude squared of the transmission, |T(w)|", as the frequency w and couplings gin and gout are varied. Here, 
Aw = cij — 0)0, where coo is the resonant frequency of each nanoparticle. Panels (a), (e) and (i) correspond to couplings gi„ = gout for an array of 
n = 3, 5 and 7 nanoparticles respectively. Panels (b), (f) and (j) correspond to gin = 2gout and panels (c), (g) and (k) correspond to gi„ = gout/2. 
Panels (d), (h) and (1) have a threshold placed on the value of |T((x))p, with the solid area corresponding to |T(w)p > 0.98 for gi„ = gout. 



Increasing (decreasing) g^p shrinks (expands) all axes. This 
observation can be useful when comparing two gnp regimes 
with each other. Note also that max(gin/g™'',^out/^™'') - 1 
must be imposed, where = O.lwo, otherwise we would 
move away from the weak coupling regime for the source 
and drain. In other words, the rescaled couplings gin/gnp and 
gout/gnp can in principle go higher than 1, but the value for 
gnp must be lower than 0.1 wo to compensate so that we are 
still in the weak coupling regime. In Fig. |4] (b), (f) and (j) 
we plot a different cross-section for n - 3, 5 and 7 nanopar- 
ticles, where = 2^out and in Fig. |4] (c), (g) and (k) for 
^in = goui/2. In Fig. |4] (d), (h) and (1) we place a thresh- 
old on the value of |T(w)P such that the solid red area corre- 
sponds to |T(w)p > 0.98 for gin = gout- One can see that as 
the source and drain couplings increase, the range over which 
the transmission is close to ideal becomes enlarged about the 
central resonance, although if the couplings are too large this 
range reduces back again. Similar behaviour can be seen for 
larger odd numbers of nanoparticles, with the central 'fork' 
area becoming narrower as n increases. These behaviors can 



be understood as follows. The early increase of gi^ enables 
the off-resonant transfer from the source to the first nanoparti- 
cle, whereas its late increase leads to strong coupling as if the 
first nanoparticle becomes the extended 'tip' of the nanotip. 
A similar argument about gout applies for the last nanoparticle 
and the drain nanotip. Thus the large gin/out implies that the 
number of nanoparticles is effectively reduced to n — 2. In 
the moderate magnitude of gin/out, we have the broad region 
of frequency a> for highly efficient transfer This observation 
will be important in our study of quantum state transfer in the 
next section. 



III. QUANTUM STATE TRANSFER 

A. Qubit transfer 

We now consider quantum information in the form of a sin- 
gle quantum bit, or qubit, transferred across a metal nanopar- 
ticle array. We write the input qubit state in the source as 
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|t/r), = a\0}^ + b\l^)^ milll], where \a\^ + \b\^ = 1, and 
|0)^ and represent the vacuum state and single plas- 
mon wavepacket in the source (at the tip), respectively. The 
wavepacket is characterized by a spectral profile ^{a>) with 
daj|^(w)p = 1. More explicitly we have 

Xoo 
daj^(cj)sl{c)\Q)^. (8) 
-oo 

Then for a given input state from the source, we take both 
the nanoparticles and drain to be initially in the vacuum state. 
Using the relation in Eq. ^ and substituting for then 
tracing out the state in the source (see Appendix C), we obtain 
the output state in the drain nanowire (at the tip) as 

Pd = {\a\^ + \b\^j dcj\^(co)\Hl ~ \najt)yO)d{0\ 

XOO 
daj^*(oj)T*{co) \0)d{U 
-OO 

Xoo 
-oo 

Xoo pOO 
doj I dw'^(w)^*(w') X 
-OO J — oo 

T(co)T*(c^')\lco)d{la>'\, (9) 

where |lw)j = li^m(w) |0)^. For perfect state transfer, i.e. 
T(oj) = e'(*-^±'^), giving |r(w)p = 1 and |7^(w)|2 = 0, the out- 
put state going into the drain nanowire, described by Eq. ([8]), 
becomes the pure state pd - \4i')d{'i''\^ where \4>')d is equiv- 
alent to Eq. dS), but with ^(w) ^\oj) = ^(w)e'(*-'±''> and 
s\j,(jj) — > dl^^{(jS). The change in the spectral amplitude of the 
wavepacket is equivalent (upon Fourier transform) to a pos- 
itive temporal shift (delay) in the wavepacket, which corre- 
sponds to the time that the wavepacket takes to move from 
the source tip to the drain tip. For concreteness, consider a 
Gaussian wavepacket with spectral amplitude profile 

^{(jj) = {27T(T^y^'''e (10) 

where ojq is the central frequency and cr - Scoj (2 v21n2) 
is the standard deviation corresponding to a FWHM band- 
width S(jL> for the spectral intensity profile |^(w)p. Applying 
the transform ^{cS) — > ^'{oS) - ^(w)e'^*'*^'^' and assuming a 
small enough so that there is linear dispersion about uo, 
then k ~ w«eff/c = <y/ceff, where n^s and Cetf are the effec- 
tive refractive index and speed across the nanoparticle array. 
We can then write kx ~ wx/Qfj- - a>6t, where 6t is the time 
taken for the wavepacket to propagate from the source tip to 
the drain tip and x is the total effective distance (from the cen- 
tre of the source tip to the centre of the drain tip). Setting 
^(w) — ?> ^'(w) = ±^(<jj)e'"* and taking the Fourier transform 
one finds 

^'(f) = ±(2cr2/;r)i/'*e~'"'^'-*''-'"''('-*' EE ^(f ~ 6t), (11) 

corresponding to a positive shift, or delay, of 6t in the time 
domain. 

We now consider the fidelity of the transfer, defined as 
F - d {^\Pd W)d^ where \4'')d is the ideal transferred state 



including the dispersion, as defined previously. The fidelity 
describes how close the output state is to the expected one, be- 
ing zero for orthogonal states and 1 for perfect transfer. Thus 
we use it to quantify the quality of state transfer. A straight- 
forward substitution gives the more explicit form 

Xoo 
du\^{ut{l~\nut + 2\T{oS)\) 
OO 

+ |Z7|4jJ doj\^{cDt\T{a>)^ . (12) 

Using the Bloch sphere coordinates a - cos(0/2) and b - 
e"^ sin(0/2) and averaging the fidelity over all possible qubit 

states F = ^ £ de f^" dcpF sinO, one finds \a\'^ 1/3, 
1^1'' — > l/3and|flp|fep — > 1/6. Thus, for a given nanoparticle 
array and input wavepacket defined by ^{oj), with a knowledge 
of |T(w)|, one can calculate the average fidelity of the output 
qubit state going into the drain nanowire using Eq. (fT2l i. Note 
that Eq. (fT2] i is irrespective of dispersion and depends only on 
|T((y)|, since we have taken the fidelity with respect to 
setting ^'(w) correctly to the expected profile resulting from 
an arbitrary input ^(w), which compensates the dispersion of 
transmission. However, for simplicity we limit our discussion 
to linear dispersion, where the expected output state by perfect 
transfer has the profile given in Eq. (fTTT l. 

In Fig. ID (a), (c) and (e) we show the average fidelity F 
for an array of n = 3, 5 and 7 nanoparticles. Here one can 
see immediately that for a small enough bandwidth, the state 
can be transferred across the array with perfect fidelity. The 
dashed lines correspond to fidelity contours, with the lowest 
curve (0.66 ) corresponding to the classical threshold for a 
quantum channel: the best fidelity achievable by measuring 
an unknown qubit along a random direction and then sending 
the result through a classical channel using classical correla- 
tions ll60ll . The solid blue curves bound the region (from be- 
low) in which the dispersion is approximately linear, so that 
we can use the approximation k ~ wneff/c = w/ceff to obtain 
the form of the expected output spectral profile ^'(w). This re- 
gion is found by calculating the group velocity vg(w), where 

'(w) - dkldio' \oj'=a} and k is found from Eq. (|7]). For linear 
dispersion about the resonant frequency we should have that 
vcito) — VG(<yo). In Fig.|5](b), (d) and (f) we show the scaled 
group velocity vg(w) = vg(w)/vg(wo) for an array of « - 3, 
5 and 7 nanoparticles. One can see that there is a wide fre- 
quency range available in the linear dispersive regime, given 
a large enough input/output coupling can be achieved. 



B. Single-photon and coherent state transfer 

We now discuss the transfer of two particular kinds of in- 
put state: single-photon states and very low-intensity classical 
light described by coherent states having an average photon 
number of 1 . These are typical quantum and classical states 
of light, respectively, and while they appear to be similar, they 
are in fact very different states altogether, with different mea- 
surable physical properties. On one hand, a single-photon 
state injected into the source nanowire can be described by 
|1^)^ = X^da;^(o.)4(6>)|0),, with ^°"^du\^iu)\^ = 1. On 
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FIG. 5: Fidelity of quantum state transfer for a single qubit wavepacket in the absence of damping and scaled group velocity over a range of 
frequencies and couplings, (a): Average fidelity F for transferring a qubit wavepacket over w = 3 nanoparticles as the bandwidth (cr) and 
coupling gin (= gout) are modified. All parameters are scaled by g„p (b): Scaled group velocity \'a{(^) = VG(<jj)/vG('jJo) as it deviates from that 
at the resonance frequency wo for n = ?>, showing the regions of approximate linear dispersion, (c): Average fidelity for transferring a qubit 
wavepacket over n = 5 nanoparticles. (d): Scaled group velocity vg(w) as it deviates from that at the resonance frequency wo for n = 5. (e): 
Average fidelity for transferring a qubit wavepacket over n = 1 nanoparticles. (f): Scaled group velocity va(io) as it deviates from that at the 
resonance frequency ojq for n = 7. In all group velocity plots, the region inside the blue lines for Aa> corresponds to the frequency range below 
the blue lines for cr shown in the average fidelity plots. 



the other hand, a coherent state is described by |{q'}) j - 
exp(l]^^ — iin,ff)|0)j, where the wavepacket operators are 

t,a = rido;ff(w).tM. with J^^daj\a(aj)\^ = («) 61. 
Using the quantum theory we have developed to describe the 
nanoparticle array system, we now highlight a difference be- 
tween single-photon states and coherent states (which are con- 
sistent with classical electromagnetic theory). The aim is to 
show the necessity of our quantum formalism in order to cor- 
rectly predict measurable physical properties of the transfer 
process. 

First we consider that the average photon number of the in- 
jected coherent state is 1, i.e., (n) = /_°^dw|Q-(w)|2 = 1, and 
the wavepacket amplitude a{cS) is the same Gaussian form as 
^(a>). The scattering matrix given in Eq. ^ enables us to treat 
the nanopartice array as an effective beam splitter, and for a 
single-photon state and coherent state we obtain the following 
respective output states at the nanotips 

^ r dJacomoj)\l^)^+^(aj)T(co)\uX 

J — OO ^ ' 

It is clear from the above that each input state arriving at the 
source nanotip is transmitted and reflected in a different way: 
the single-photon state becomes an entangled state of trans- 
mitted and reflected single-plasmon states while the coherent 



state remains as a separable state of transmitted and reflected 
coherent states of plasmons. Nevertheless, the detection prob- 
abilities (mean excitation flux) at the drain are exactly the 
same as each other. This is calculated by finding the expec- 
tation value («£/„,„), where n^^^, - da)dl^^((L))dout{<^), and 
gives the same result for both input states 



(«rf„J= dc^|^(w)|2|r(w)|2 = dcj\a(aj)\^\ncj)\\ 

This implies that there is no difl^erence in the energy trans- 
fer efficiency between single-photon states and coherent states 
when they are injected into the nanoparticle array. However, 
in quantum information processing, and in particular quantum 
communication, a more meaningful measure of the transfer 
success is not the energy efficiency, but how well the informa- 
tion content that is encoded into a physical state is preserved. 
This can be quantified by the fidelity between the transfeiTed 
state and the ideal transfeiTed state, as defined in the previous 
section and it is a measurable physical property of the transfer 
process; it can be measured by performing quantum state to- 
mography 1 63]. The fidelity for the transfer of a single -photon 
state is obtained by substituting a - Q and ^7 = 1 in Eq. ( fT2b . 
The fidelity for the continuous-mode coherent state transfer 
is summarized in Appendix D. The respective fidelities are as 
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follows. 



exp 



XC30 
daj\^(co)\^\n 
OO 

Xoo 
dc,\a{ujt(\n 
-OO 



and 



It is clear that they are not the same. It is important to note that 
while the transfer of a single-photon state and coherent state 
are equivalent in the sense that the nanoparticle array trans- 
mits the same amount of their energy from the source to the 
drain nanowire, they are in fact different from the viewpoint 
of the transfer of information encoded within the states. This 
behaviour naturally carries over to the general case of qubits, 
where a and b are arbitrary, as in Eq. ( fT2l l. It also applies when 
damping is introduced (see next section). 



IV. DAMPING 

A. Physical model and transmission 

We now include damping in our model. The effects of loss 
in the system are due to the interaction of the electrons (sup- 
porting the surface plasmon field) with phonons, lattice de- 
fects and impurities |34, 61], as well as radiative scattering of 
the surface plasmon into the far- field [34]. For most scenarios 
of nanoparticle arrays, the couplings between nanoparticles 
are large enough such that most of the field remains within 
the array, with radiative scattering rates generally 5 orders of 
magnitude smaller than the relaxation rate |34]. Thus we as- 
sume radiative scattering can be neglected in our model. This 
assumption also allows us to neglect possible scattering at the 
tips. Electronic relaxation effects on the other hand cannot 
be neglected and lead to damping of the supported surface 
plasmon field. In our model we describe this as an amplitude 
damping channel at each nanoparticle. In this context a mech- 
anism can be introduced where the damping is modeled by 
coupling of the field at each nanoparticle to an independent 
bath mode, which is eventually traced out from the system 
dynamics, as shown in Fig. |6] As we are interested in the 
mapping of the input field at the source tip to the output field 
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FIG. 6: The plasmonic nanoparticle array including bath modes to 
model damping at each nanoparticle. 
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FIG. 7: Transmission spectral profiles and dispersion relations for 
damped arrays of nanoparticles. Panels (a)-(e) correspond to the 
amplitude squared of the transmission, |T((jj)|", as the frequency is 
varied for an array of m = 1, 2, 3, 5 and 7 nanoparticles respec- 
tively. Panels (f)-(j) correspond to plots of the real part of the ef- 
fective wavenumber k' for n = 1, 2, 3, 5 and 7 respectively. Also 
shown are points corresponding to the k''^ transmission resonance 
peaks from panels (a)-(e) as well as the dispersion relation for the in- 
finite array case (dashed line). Here the couplings used are the same 
as the undamped case, i.e., g„p = — O.lwo and gin/out = 0.01 wo. Note 
that larger transmission values can be achieved by increasing these 
couplings, as explained in the text and shown in Fig. [8] 



at the drain tip, we assume that the source and drain excita- 
tions experience no loss when propagating in/out of the tip 
regions. Such insertion loss can however be incorporated us- 
ing standard waveguide methods (SIS^ although a specific 
model will depend on how the fields in the nanowires are ex- 
cited and collected, for instance, how far they propagate in the 
nanowires. Various types of dielectric-metal structures can 
significantly reduce these losses ll62ll . 



10 



(ajn = 3 {b)n = 5 (c)n = 7 „, ,2 




■■272 4V2 " Ve 2V6 ° 4/V3 a/V3 

9inl9np = 9oui/9np ffin/ffnp = Souf'fnp Sin/ffnp =5ojt/ffrip 



FIG. 8: Damping for the amplitude squared of the transmission, |T(oj)|^, as the frequency oj and couplings gin and gout are varied. Panels 
(a), (b) and (c) correspond to couplings gi„ = gout for an array of n = 3, 5 and 7 nanoparticles respectively. 



The scattering matrix in the presence of damping is derived 
in Appendix E. In the forward direction, we have the relation 
between the input field operator sl^icij), the output field opera- 
tors s Jut^*^) ^^'-^ '^out^'^)' '^^'^ ^^'■h operators A^^^^ .(w), 

n 

i=l 

(13) 

where the index s is dropped in the coefficients for ease of no- 
tation. The i-th nanoparticle loss coefficients, Siiai), are also 
functions of the system parameters ^in, gout, gij and lo,, and 
|7^(w)p + |r(w)p -I- E"=i \Si(to)? = 1- This method again 
allows us to describe the nanoparticle array as an effective 
waveguide, with T(w) = |T'(w)|e'^'^-*^'^', where |T(w)| is the 
transmission in the ideal case (no damping) and the wavenum- 
ber k = k'' + ik' has become complex as a result of the damp- 
ing [I44I], which now depends on the system parameters g\n, 
gout, gij, and the relaxation rates F,- at each nanoparticle. 
Thus, we have that \T (w)| = \ f (u)\e-'''\ 

In Fig.|2](a)-(e) we plot the amplitude squared of the trans- 
mission, |T(<jj)p, as the frequency is varied for an array of 
« = 1, 2, 3, 5 and 7 nanoparticles respectively. To compare 
the damping with the ideal case shown in Fig. [3] we use the 
same system parameters: all local frequencies are equal - 
Wo, V/, the couplings are equal gi j = g^p = — O.Iwq, V/ and its 
nearest neighbors /', and the source and drain couplings are set 
as gin/out = O.OIwq. In Appendix F we provide the analytical 
form for the T(<y)'s with damping. The damping rate for each 
nanoparticle depends on its size and is given by Matthiessen's 
rule II34II : F = vf/zIbH-v/t/^, where for silver /1b = 57nmisthe 
bulk mean-free path of an electron, Vf = 1.38 x lO^m/s is the 
velocity at the Fermi surface, and the effective radius R ^ R. 
We use Fj - 0.0158wo,Vi, which corresponds approximately 
to the damping rate for a silver nanoparticle with a radius R in 
the range 20— lOOnm. For a given «, the transmission spectral 
profiles in Fig. |7](a)-(e) again have n resonances at frequen- 
cies (L),y - <jj() + Igcos.if'^d), where k''^ - jn/{n + l)d for 
j = However, the width of the resonances has been 

broadened and the height lowered as a result of the damping. 
In Fig.|2](f)-(j) we plot the real part of the effective wavenum- 
ber k'' from Eq. (Q for « = 1, 2, 3, 5 and 7 respectively. 
Note that Eq. O remains valid, as the imaginary part of the 



wavenumber k is absorbed into the magnitude of the trans- 
mission, |T(w)|. Points corresponding to the k'''-' transmission 
resonance peaks from Fig.|2](a)-(e) are marked. Also included 
in these figures, as before, is the dispersion relation for the in- 
finite array case (dashed line). 

From Fig. [T] (a)-(e) one can clearly see that the transmis- 
sion peaks are much reduced from the ideal values. How- 
ever, despite this, it is possible to increase the maximum peak 
value by increasing the source and drain couplings, as shown 
in Fig. [8] (a), (b) and (c), where we plot a cross-section of 
the amplitude squared of the transmission, |T(w)p, as the fre- 
quency a> and coupling gin {- gout) are varied for an array of 
« = 3, 5 and 7 nanoparticles respectively. Here, Aw - a> — a>() 
and as before, all parameters are in units of the nanoparticle 
coupling gij = gnp, V/ and its nearest neighbors j. One can 
see from Fig. [8] that as the source and drain couplings (gin and 
gout) increase, the transmission maximum can be increased, 
although ultimately the damping dominates the transmission 
as n increases, as can be seen by comparing Fig.[8](a) with (c). 



B. Qubit transfer 

We now discuss the fidelity of state transfer for a single 
qubit wavepacket state under realistic conditions of loss at 
each of the nanoparticles. After including the bath modes at 
each of the nanoparticles, one finds that the expression for the 
fidelity given in Eq. fT% remains valid (see Appendix C), with 
the fidelity depending only on the absolute value of transmis- 
sion coefficient. In Fig.|9](a), (c) and (e) we show the aver- 
age fidelity F for an array of « = 3, 5 and 7 nanoparticles. 
The dashed lines correspond to fidelity contours with the low- 
est curve (0.66 ) corresponding to the classical threshold for 
a quantum channel, as before. The solid blue curves bound a 
region (from below) in which the dispersion is approximately 
linear, vg(w) ~ vg(wo). In Fig.|9](b), (d) and (f) we show the 
corresponding scaled group velocity vcia>) = vdcS) / vdcjo) 
for n = 3, 5 and 7 nanoparticles. For n - 3, one can see 
in Fig. |9] (a) that the nanoparticle array can provide a trans- 
fer channel giving an average fidelity of up to ^ 0.93 even 
when damping is present, in which for large bandwidths cr the 
source and drain couplings gin and gout must be increased to 
values close to the limit of the weak coupling approximation. 
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FIG. 9: Fidelity of quantum state transfer for a single qubit wavepacket under damping and scaled group velocity over a range of frequencies 
and couplings, (a): Average fidelity F for transferring a qubit wavepacket over n = 3 nanoparticles as the bandwidth (cr) and in/out couplings 
gin = gout are modified. All parameters are scaled by g„p (b): Scaled group velocity vo(w) = VG((jj)/vG((jJo) as it deviates from that at 
the resonance frequency aio for n = 3, showing the regions of approximate linear dispersion, (c): Average fidelity for transferring a qubit 
wavepacket over n = 5 nanoparticles. (d): Scaled group velocity vaiijj) as it deviates from that at the resonance frequency ojq for n = 5. (e): 
Average fidelity for transferring a qubit wavepacket over n = 1 nanoparticles. (f): Scaled group velocity vg(w) as it deviates from that at the 
resonance frequency uiq for n = 7. In all group velocity plots, the region inside the blue lines for Aw corresponds to the frequency range below 
the blue lines for cr shown in the average fidelity plots. 



I^in.out/^™"! = 1. Note that in these plots one cannot decrease 
the coupling gnp in order to reach gin,out values much larger 
than 1, as we have set gnp = — O.lwo, unlike the ideal case 
where it could be modified. The reason for this restriction is 
that reducing the nanoparticle coupling gnp means the damp- 
ing rates begin to dominate, lowering the maximum transmis- 
sion and average fidelities further as a result. 

For n - 5, one can see in Fig. |9](c) that the maximum av- 
erage fidelity attainable is ^ 0.88; no contour can be plotted 
for 0.9 or above, regardless of the bandwidth cr. For n - 1 
and above, this situation then becomes gradually worse and 
one can see in Fig. |9](e) that although the maximum average 
fidelity attainable is ^ 0.84, the source and drain couplings 
need to be increased close to the weak coupling limit, in addi- 
tion to the use of a narrow enough bandwidth. 

The results obtained here indicate that only small-sized ar- 
rays with « < 7 are useful for the transmission of qubit 
states encoded into the number state degree of freedom. How- 
ever, it may be the case that for particular applications, short- 
distance communication (< /vm) is required at optical frequen- 
cies, making the use of a nanoparticle aiTay quite beneficial. 
For example, the nanoparticle waveguide could be used as an 
enhanced mediator between emitter systems on a very small 
scale. On the other hand, additional degrees of freedom for 
the LSP excitations, the embedding of emitter systems into 
the waveguides, novel types of metals with reduced damping 



rates and new schemes for achieving gain in plasmonic me- 
dia may enable one to eventually counter the effects of loss 
highlighted here. 



C. Plasmon interference 

In Section III B we showed that in order to correctly de- 
scribe the transfer of a quantum state through a metal nanopar- 
ticle aiTay one requires the quantum formalism we have de- 
veloped in this paper. Here, as an additional example of the 
necessity of a quantum formalism for the metal nanoparticle 
array, we investigate the interference of two plasmons. We 
consider the plasmons enter the array from opposite ends, one 
from the source and the other from the drain nanowire. The 
input state at the nanotips in this case can be written as 

Xoo poo 
dw, I Aud^{(i),,u)d)s\j,u,;)d\^i(Dd)\0) 

where |0)^^^ denotes the vacuum state for the source, 
drain and baths and the normalization of the state 
vector imposes a nomalization on (/'(w,, w^/), so that 
^AlOs j_^A(ji)d\il/{ci>s,^d)\ - 1- By using the scattering 
matrix given in Eq. (IE-2b of Appendix E, that describes for- 
ward and backward propagation in the aiTay, we have the out- 
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FIG. 10: (a): Plasmon interference. Here plasmons enter the nanoparticle array from both sides. One plasmon from the source and another 
from the drain. Nonlinear absorption occurs via quantum interference, even though the damping in the array is linear, (b): Survival probabihties 
of zero-plasmon (yellow), one-plasmon (red), and two-plasmon (blue) for n = 3, as the coupling g;,, (= gout) is varied for Aa> = 0. (c): Survival 
probabilities for n = 3 when the one-plasmon survival probability is minimized as the loss T is varied for Ato = 0. The values of gin (= gout) at 
which minimization occurs are shown in the inset, (d) and (e): Reflection and transmission coefficients corresponding to the couplings in (b) 
and (c) respectively. 



put 

l'A)out = 

(7^,(«,)4,(w,) + TsiojMljco,) + Fl{co,)j X 

where the noise operators are defined as Fl(a>s) = 
2,-5j,,(Ws)A^,.(Ws) and ^^(wrf) = Y,iSd,iiujMl^i/(^d)- We 
consider small bandwidths for i/'(ws, ojd) over which the trans- 
mission, reflection and damping coefficients do not vary ap- 
preciably and therefore these coefficients will be approxi- 
mated as frequency independent. The case of gin - gout 
is considered, so that we have 7^(w) = Ts(<y) = T and 
TZdico) - Ti-Apj) - Tl- The probability of finding two plas- 
mons in the source nanowire and none in the drain nanowire 
is then (see Appendix G) 

p(2,,0rf) = |7e|2|rp(i+i), (14) 

where we have introduced the (real) overlap integral 

XC30 poo 
doJs I d(jJd4'{(^s,(JJd)4'*i.^d,OJs)- (15) 
-OO — OO 



Here, unit quantum efficiency of the photon detector and in- 
finite counting time are assumed. Similarly, the remaining 
nonzero probabilities are 

P(0,,2d) = P(2.,,0d), 

p(i,,irf) = \n\U\T\U(n^T*^ + n*^T^)i, 
p(i,,Orf) = (|7^|2 + |rp)(i - 17^|2 - in') 

-iUT* +'R*TfX, 

P(Os,U) = P{ls,Qd\ 

p(o,,Orf) = {i-\n\^ ~\T\^f + {n*T + T*n)-T. 

(16) 

Here, for simplicity, we consider that the plasmons have the 
same wavepacket profile, i.e., iff{ciJs,a>d) - ^itiJs)^ia>d), where 
^{a>) is given in Eq. ( fTOl i. If I = 0, the Fourier transform 
of the spectral amplitudes for the two plasmons do not over- 
lap in time in the nanoparticle array, and the probabilities in 
Eqs. (fT4b and ( fT6l ) describe the case of two independent parti- 
cles 1 64]. On the other hand, if I = 1, the Fourier transform 
of the amplitudes overlap perfectly in time. In this case, tem- 
poral and spectral indistinguishabilities are immediately satis- 
fied and for \TZ\^ - |Tp = 1/2 one recovers the well-known 
Hong-Ou-Mandel (HOM) quantum interference effect l65ll . 
where the two excitations are always found to be in the same 
output mode: P(2„ Qd) = P(Oj, 2d) = 1 /2, with all others be- 
ing zero. In general, however, when \TZ\^ + |Tp or damping 
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is present (\Tl\^ + \T\' < 1), the probabilities for two, one, 
or no plasmons to survive are f 2 = ^(2j, 0^) + P{Qs, 2^) + 
P(l,, Id), Pi = P(h,Oa) + P{Os, I A and Pq = PM^), re- 
spectively, with Y^i^qPi - 1. 



For an array of n = 3 nanoparticles, we plot in Fig.[TOl(b) 
the survival probabilities as the coupling gin (= gout) is var- 
ied for Aw = 0, r = 0.0158(^0, and gnp - — O.lwo. One can 
see that the probability for one of the plasmons to survive (or 
be absorbed). Pi, can be very low depending on the in/out 
coupling. Indeed, at a particular point marked by the dashed 
line, nonlinear absorption occurs: either both plasmons are 
absorbed, Pq ^ 1/2, or neither is absorbed, P2 ^ 1/2, 
even though the damping in the nanoparticle array is a lin- 
ear process. Surprisingly there is no one-plasmon absorption. 
Pi ^ 0. This effect is due to quantum interference of the plas- 
mons and cannot be described in terms of a classical treatment 
of the nanoparticle array ll64ll . In Fig. [IO](d) we show the 
corresponding reflection, |7?,p, and transmission, |Tp, coef- 
ficients. The transmission coeflicient in this plot can also be 
seen by taking a cross-section from Fig.|8](a) at Aw = 0. One 
can see in Fig. [10] (d) that the nonlinear absorption effect is 
maximized at a similar point to that for the HOM interfer- 
ence effect: reflection and transmission coefficients are equal- 
ized, but at 1/4 instead of 1/2 due to the necessary presence 
of damping in order to see nonlinear absorption ll64ll . 

In Fig. [To] (c), we show how increasing the loss at each 
nanoparticle affects the two-plasmon interference for « = 3. 
Here, Pi is minimized by modifying gia {- goad as the loss F 
is increased for Aw = 0. The corresponding Pq and P2 are also 
shown. One can see that nonlinear absorption can be made to 
occur over a large range of loss. The values of gi,, (= gout) at 
which Pi is minimized are shown in the inset and the corre- 
sponding reflection and transmission coefficients are shown in 
Fig. [13(e). Note that as the amount of loss increases, both the 
minimum value of Pi and the required coupling gin (= gout) 
are increased also. In particular, one can see in Fig. [TOl (e), 
that as the damping in the array increases, it becomes impos- 
sible to equalize the reflection and transmission coefficients 
by changing gin, as the transmission is affected more by loss 
within the array. This asymmetry leads to an eventual break- 
down of the quantum interference effect and subsequently the 
nonlinear absorption. 

The behaviour shown in Fig.[TO](c) allows us to predict the 
growing trend of the minimum value of Pi and the optimal 
value of gin (= gout) as n increases. This is because the overall 
amount of loss in the array effectively increases as the number 
of nanoparticles is increased. For Aw = 0, F = 0.0158wo and 
gnp = -0.1 wo, one finds that Pf" = 0.012, 0.034 and 0.063 
when g°J'Vgnp = 0.1543, 0.2223 and 0.2824 for n = 3, 5 
and 7, respectively. The corresponding zero and two-plasmon 
probabilities are Pq = 0.4999, 0.4995 and 0.4990, and P2 = 
0.4880, 0.4663 and 0.4380. Thus, nonhnear absorption by 
two-plasmon interference is present in the nanoparticle array 
for M = 3, 5 and 7. The nanoparticle array may therefore act as 
an effective two-plasmon absorber, despite the linear optical 
properties assumed in the model. 



V. SUMMARY 

In this work we studied the use of an array of metallic 
nanoparticles as a channel for on-chip nanophotonic quantum 
communication. After introducing the model for the physical 
system in the quantum regime, the transfer of a quantum state 
encoded in the form of a single-qubit wavepacket was studied 
under ideal conditions. We then showed the necessity for our 
quantum formalism in predicting the outcomes of measurable 
physical observables. The effects of loss in the metal were in- 
cluded in our study, thus putting the investigation into a more 
practical setting and allowing the quantification of the perfor- 
mance of realistic nanoparticle arrays as quantum channels. 
For this task we used the average fidelity for the state trans- 
fer We found that small-sized arrays are practically useful for 
the transmission of qubit states encoded into the number state 
degree of freedom. We also showed that nonlinear absorption 
can occur by quantum interference, where two plasmons are 
absorbed or neither is absorbed. Thus, the nanoparticle array 
can act as an effective two-plasmon absorber, and the obser- 
vation of this quantum interference effect may open up new 
kinds of plasmonic interference experiments in the quantum 
domain. Our study highlights the benefits as well as the draw- 
backs associated with nanophotonic periodic quantum sys- 
tems that use surface plasmons. The techniques introduced 
in this work may assist in the further theoretical and experi- 
mental study of plasmonic nanostructures for quantum control 
applications and probing nanoscale optical phenomena. 
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APPENDIX A 

Here we use input-output formalism ll52i Issll for the 
nanoparticle array to obtain an effective scattering matrix. We 
start with the Heisenberg equation of motion for an operator 
O, given by ^ = -~j,[0,H], and substitute the Hamiltonian 
H in Eq. ([T]| to obtain the equations of motion for each of the 
system operators 



dsioj) 

dt 
da I 
dt 

da, 
dt 
da„ 
dt 



-/wi(w) +gi„(w)fli, (A-1) 
--[ai,H„p]-j t/wgin(w)i(w), (A-2) 



— = ~-[ai,H„p], i-2,...,n-\. 



(A-3) 



da„ _ i . r 
— - --[a„,Hnp]- 



t/wgout(w)(i(w), (A-4) 
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dd{(jj) 
dt 



- ~iu)d{co) + ^out(w)fl„. 



(A-5) 



Here we have introduced an explicit time dependence in the 
frequency space operators s(io) and d{u)) in the source and 
drain respectively. This is because the internal field of the 
nanoparticle array may acquire some non-trivial dynamics 
which forces the external fields in the source and drain to 
have a time dependence that is different from the free field 
dynamics |52, 53]. With the above set of coupled equations 
of motion we find boundary conditions for the system before 
proceeding to solve them. Using the following solutions for 
the first and last equations (Eqs. ( I A- II ) and ( IA-51 )) 



where fo < t, with SQicS) and dQ{(jS) as the operators for s(w) 
and d(u)) respectively at time f = fo as initial boundary condi- 
tions, one finds the equations of motion (Eqs. (IA-21) and (IA-4b ) 
for the first and last nanoparticle become 



da I 

dt 

da„ 

dt 



n 
n 



H„ 



H„ 




where we have defined the input field operators as 
Sin(f) = ~(27t)~^'^ J^^ddje-'^^'-'o^soicj) and Ji„(f) = 
— (2;r)~'^^ da>e^"^^'~'"^di)(oj). Here we have assumed the 
couphngs gin(w) and gout(w) are constant over a band of fre- 
quencies about the characteristic excitation frequency being 
considered, = gin/27T and gl^^ito) = gout/^n. This as- 

sumption is valid for negligible change in the similarity of 
the modefunction profiles at the tip and nanoparticles over the 
bandwidth. We assume this can be achieved given a narrow 
enough band of frequencies along with an optimized nanotip 
geometry. 

Using alternative solutions for the first and last equations of 
motion 



s(oj) = e-'"^'-'''si{co)- ginioj) J e-'"^'-''^ai(t')dt', 
d(cj) = e-'"('^''yi(w)-§out(w) r e-'"^'-''>a„(t')dt', 



where ti > t, with si(aj) and di(cLi) as the operators for s(u) 
and d{a>) respectively at time f = f i as final boundary condi- 
tions, one finds the equations of motion (Eqs. (IA-21) and (IA-4b ) 
for the first and last nanoparticle become 



da I 

dt 

da„ 

dt 



r" tr 1 , 6in ^ 

'-[ax,Hnp\ + —ax 
n 1 

' r" i> 1 , Soul „ 
'-[an,M„p\ + —2^^n 



(A-8) 
(A-9) 



where we have defined the output field operators as 
soutit) = (2;r)-i/2j- and d^t) = 



Taking Eq. (lA-SI l and subtracting Eq. ( IA-6I 1 gives the 
boundary condition 



ai(f) = — ^(iin(f) + Sout(f))- 



(A-10) 



Similarly, taking Eq. ( |A-9l l and subtracting Eq. ( IA-7b gives the 
boundary condition 



a„(t) 



1 



-JdUt) + doutit)). 



(A-11) 



Note that throughout we assume the dispersion is negligible 
for the initial/final excitations of the source and drain fields. 
In this sense we are interested only in the relation between 
the input/output propagating fields in the nanowires near the 
tips. The dispersion during propagation of the excitations in 
the nanowires can be incorporated into the model by using 
standard methods |44]. On the other hand, the dispersion in 
the array is included in the model automatically, although we 
will need to ensure later that minimal broadening of the band- 
width due to dispersion occurs during the propagation for the 
relation gl^^fiat) = goutl^-n to still hold. We will see that this is 
a reasonable assumption for small-sized arrays. 

We use the relations [a\,H„p] - fiajiai + %i,2fl2, 
[an,H„p] - h(Aj„a„ + ngn^i,nan-i and [ai,Hnp] - fttOiOi + 
figi-\jOi-i + figi,i+\ai+\ (for / = 2, . . . , n — 1) as well as defin- 
ing the Fourier components of the nanoparticle field operators 
as fl,(f) - (27r)~'^^ e~"^'ai(cij)d(jj, V/, and we rewrite the 

source/drain operators as Sin(f) - (27r)~'^^ e'"^' Si„(cj)daj 

and din(t) - (2n)~^^^ e~"^'din(cij)dcL) (sin(a>) and din(oj) 
are general spectral operators). Then, we find upon substi- 
tution into Eqs. ( IA-3b . ( IA-61 ) and ( IA-7b the following set of 
coupled equations for the frequency operators 



i{a> ~ a>i) 



fli(w) = igiaaiico)- \/gi^Sm{tj), 



i(0J — LL>n) 



iito - tOi)ai(io) = igi-ijOi^iito) + igMjai+iito), 

for / = 2, . . . , n — 1 , 

^out' 



a„(aj) = /g„_i,„fl„_i(dL») - ^/^(iin(w). 



as well as boundary conditions from Eqs. ( lA-lOb and ( lA-llb 



:(iin(w) + iout(<y)), 



-.(dm(u) + liout(w)). 



fll(w) 

a„(cj) 



Using the above set of coupled equations we can eliminate the 
internal nanoparticle operators a, lf52[ issil to obtain 

iin(w) = 7^*(w)^out(w) + 7;*(w)Jout(w), 
dinioj) = 7^*(w)iout(w) + 7^;;(w)^iout(w), 

where the transmission Ts.d and reflection TZs,d coefficients are 
functions of the system parameters gin, gout, gij and w,-, and 
the relation |7^s,rf(w)p + \Ts.d{<^)\^ = 1 holds. 
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APPENDIX B 

Here we provide the analytical forms for the T((y)'s for 
n = 1, 2, 3, 5 and 7 nanoparticles respectively (no damping), 
where we have set ^out - gin^ 



Ti(oj) = 



i{u) — (jJq) 



4gm + (^in - 2/(w - W()))2 ' 



J--i(a)) = ~^gl„gm[{gm - 2;(w - wo))(4^^p 



-(W - W())(/gin + 2{(jJ - wo)))] 



TiCw) = 4^np5in[(2gnn(§in - 3/(w - Wo)) 



-{gin - 2i{(L) - wo))(w - Wo) ) X 
77(w) = -4g^pgin[(gnpfen - 4/(w - Wo)) 



(2,gnp - (^^ - Wo)(/gin + 2(W - Wo)))] 



-(gin - 2/(w - Wo))(w - Wo) ) X 

{Agl^ + (W - Wo)^(/gin + 2(w + wo)) 

- g^p(w - wo)(3/gi„ + 8(w + Wo)))]"'. 



APPENDIX C 

Here we show how to obtain the output density matrix for 
the qubit state entering the drain nanowire. This is done in the 
general case of damping (see Section IV). To obtain the case 
of no loss, simply set 5,(w) = 0, V/. Starting with the single 
qubit wavepacket in the input modes of the source nanowire 



dw^(w)^|'„(w) |0) , 

-oo 



and making use of Eq. (fT3T l of Section IV, (equivalent to 
Eq. (|5]l of Section 11, when iS,(w) = 0, V/), substituting for 
s\j.oj) one obtains the state Ps,d,A - \'P) s,d,A{'P\ which describes 
the total state in the external output modes of the source- 
nanoparticle-drain system, where 

Xoo p 
dw^w) 7^(w)4.(w) 
-oo L 
n 

+r(w)JL(w) + Si(aj)Al,/c,) |0) ,_, 



/■=! 



Removing the source modes from the description of the state 
Ps,dA is achieved mathematically by tracing them out to give 

Xoo 
dw j(l„|p.v,rf,A lU).;- 
-OO 

Tracing out the / bath modes recursively in a similar way gives 
Pd = (\a\^ + \b\^ j dw|^(w)p(l - |r(w)|2)) |0)X0| 

XOO 
dwr(w)r*(w) \o)d{U 
-oo 

Xoo 
dwaw)r(w)|i„)xo| 
-oo 



+a 



+ \b\ 



f dw r dw'^(w)f*(w') X 

— OO "J 

noj)T*(cj')\i^)d{K'\. (c-i) 



APPENDIX D 

Here we derive the fidelity for the coherent state transfer 
described in Section 111 B. First, we consider two continuous- 
mode coherent states defined as 

\{J3}) = exp(^;-V|0), 

where the photon wavepacket operators are given by = 
dwa(w)^^(w) and - dwyS(w)£^^(w), with 

J^dw|Q;(w)p - ha and J^dw|yS(w)p = hp. The opera- 
tors ^^(w) {b{<jj)) represent the creation (annihilation) opera- 
tors associated with a field excitation, which obey the bosonic 
commutation relation {b{<jj),b^{(Jy\ - 6{(jd — cJ). In general, 
the fidelity between the two continuous-mode coherent states, 
defined as f = | ({yS}|{a}) p, is found by direct substitution 
to be 



F - exp[- 



f 



dw|a(w)-je(w)p]. 



The fidelity between the transferred coherent state HaT}) 
in Section 111 B and the ideal transferred state lla'}), where 
a'(a>) - Q'(w)e'*'*, is then 



F — exp 



f 



dw|a(w)|2(|r(w)| - If 



APPENDIX E 

Here we use input-output formalism for the nanoparticle 
array under realistic conditions of loss at each nanoparticle, 
where the damping is modeled by coupling of the field at 
each nanoparticle to an independent bath mode. The inter- 
action of each nanoparticle to its bath mode takes the same 
form as the coupling of the first nanoparticle to the source 
nanowire tip, except with a coupling strength determined by 
the rate of damping to match the classical case, as done pre- 
viously for the inter-particle couplings gi j. This approach as- 
sumes a weak d amping rate and Markov approximation for 
the bath modes 1521 15311 . To mathematically incorporate the 
bath modes into our model, the original coupled equations are 
modified to become 



!(w-wi)- Y - Y 



/(W - W,) - — 



ai(w) = /gi,2fl2(w) - ^/ik^in(w) 

- VriAin,i(w), 

a,(w) = /g,_i,,a,_i(w) + igi+iiCti+iiu) 

r,Ain,;(w), for i - 2, . . . ,n — 1, 



16 



1{(AJ - W„) y 



a„(aj) = /g„_i,„a„_i((y) 



where F, corresponds to the electronic relaxation rate at 
nanoparticle /. Extra boundary conditions are then imposed 
on the system dynamics given by 



1 



=(Ai„,,(w) + Aouu(w)), for ! = 1, . . . ,« (E-1) 



As in the case of the source and drain nanowire system, the 
bath operators A{a>) obey the bosonic commutation relations 
[A{a>),A^{a)')] = 6{u) — <J). By solving the above new set 
of coupled equations and eliminating the internal nanoparticle 
operators a,, we obtain a scattering matrix linking the source, 
drain and bath operators. Then for a given input state from 
the source, we take the bath modes (and drain) to be initially in 
the vacuum state |0). The scattering matrix is applied and the 
bath modes are traced out to obtain the effective transmission 
in the array. In general one finds 

slico) = 7^,(w)4^(«) + Ts{oj)dljpj) + Yj Ss,i{u))Al^,j{0J\ 

i 

(E-2) 

where |7e,,^(«)|2 + \Z,d{c^)\^ + 2/ = 1- 

APPENDIX F 

Here we provide the analytical forms for the T((y)'s for 
n = 1, 2, 3, 5 and 7 nanoparticles respectively with damp- 
ing, where we have set ^out = ^in. 



TIM 



2gm + r - 2i{(i) - wo) 

in 

4^np + (^in + r - 2;(w - wo))2 ' 

-8gnp^in[fen + T - 2/(w - Wq)) X 

igl^ + (T~ 2i{u - wo))fe„ + r - 2/(w - wo)))]"', 



TsM = 32g^^gin[(4g2^ + (r + 2/(w - wo)) X 

{gin + r - 2/(w - Wo))) X 

(Agl^{2g,, + 3(r - 2/(w - wo))) 
+ (r - 2/(w - wo))2(gi„ + r - 2/(w - Wo)))]- 



and 



n^gl^iAmUgin + 2(r - 2/(w - wo))) + 



5np 



(r - 2!(w - wo))^fen + r - 2/(w - Wo))) X 
(32gl + 4glJ3gi, + 4(r - 2/(w - wo))) x 



(F - 2/(w - wo)) + (F - 2/(w - wo))-^ x 
fe„ + F-2/(w-wo))))]-'. 



APPENDIX G 

In our discussion of two-plasmon interference we needed 
to evaluate the probabilities for finding plasmons in the output 
state given in Eq. (fl4l i. For an arbitrary state \tfr), the prob- 
ability of detecting one plasmon in a given mode at any fre- 
quency w is given by P(l) - J_ dw| p, where |n„) 
is the continuous mode number state, as previously defined. 
The probability of detecting two plasmons in a given mode, 
one at any frequency w and the other at any frequency w' is 
thenP(2) = X°^dwX°^dw'|(2(^,<^ol^^)l^ where |2(„,„,)) = 
■^'^out('^)'^out('^') |0) is the continuous mode pair-state jiill . 
which allows for each plasmon to have a different frequency 
profile. Thus we have the following probabilities 



P{2s,Qd) = 

Pi0s,2d) = 

P(ls,Od) = 

PiOsAd) = 

PiOs,Od) = 



- I da> I dco' X 

— OO "J — oo 

I s,dA (0| ■Sout('y')iout(w) \tf))^^i P, 

- I doj I da>' X 

J — oo J — oo 

I s,dA (0| domi(^')domi(^) |i/')out 1^ 

Xoo poo 
doj I doj' X 
oo "J — oo 

I s.dA (0| ■Sout(^^)'?out('^') |lA)out 1^ 

Xoo 
dw| (0| .Sout(w) |lA)out l^ 

-OO 

Xoo 
dw| (0|^iou,M|<A)ou,|', 
-oo 

I s.d,A (0|«A>out I'- 



These may be evaluated by using the relationship between the 
input and output field operators given in Eq. (IE-2l i and the 
following commutation relations 



[iout(w), 

[iout(w), ^?o^^t('^)] 
[iout(w),/'J(w)] 



6((xJ ~ w') = [Jout(<y), '^L('^)] ^ 
[(iout(t^), Jout('^)] = 0' 
[.?out(w), f\{oS^ 
[iout(w), A(w)] 
[iout(w), FrfCw)] = 0, 



and 



[f,(w),fJ(w')] 
[/■^(w),/-](w')] 

[A(t^),f'l(w')] 

[f^(w),fJ(w')] 



5(w 
5(w 



w')(i - |7^,(w)p - ir.(w)p), 
w')(i - |7^^(w)|2 - |r/w)|2), 

— 5(w — w')(7?.*(w)7rf(w) -I- T^{u)TZd{oS)), 

-s{u - w')(r;(w)7^,(w) + 7^*(w)7;(w)). 



By straightforward substitution, one finds 
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*J — oo *J — oo 

P(2s,0d) = - I I dwXw,w')7;M7ed(a>') + ^i'(w',a>)7;(a>')7ed(w)p, 

*J — OO ^ — oo 

XOO /^OO 
dw dw'|iA(w, oj')ns((^)nA(o') + ^(pJ, oj)Ts{aj')Td{cj)\\ 
- DC' J — OC' 

duo dw'maj,(o')\^\nAcj)\\l - \nd(<o')\^ - \Td((o')\^) 

■oo *J — oo ^ 

dw dw' |iA(w,w')|2|rv(w)|2(l - \Tld{Jt - XTdioJt) 

-CX) *J — oo ^ 

- i//*(co,u')ifr(j,c^)T;{(^)nd(co)(T;(co')nAco') + n*j(J)Z(co')) 
- riJ,ojmco,J)n';iioj)rsioj)CJz;ioj')rdioj') + r;ioj')ndioj'))j, 

XOO /^OO - 

dco da)'ma>,co')\\l - |7e,M|2 - \Ts(o>)\^)(l - \Tld(co')\^ - \Td(o>'t) 
■OO *J — oo ^ 

+ r(o^', ojMoj, J)iTd*(oj')n,(oj') + ni{J)Ts{oj')){nl{J)Td{oj') + T:{aj')nd{J))\, 

I 



where only tenns making a nonzero contribution have been P(0,, \d) 
retained. Considering the case where the transmission coeffi- 
cients are approximately constant over the range of frequen- 
cies for which \il/{a), oj!)\\s significant one finds 



P(2„0d) « \ns\^\Td\\l+T), 
Pi0s,2d) « \Ts\^\nd\\l+T,, 
PihAd) « in/iUdl' + \%\^\Td\' 

+{n,ndT;T; +n:n*Tsrd)x, 
Pih,Od) « \n,\\i - lUdl^ - \Td\^) 

+\rd\\i-\nsf-\rsf) 

-(2|7^/|rd|' + KKTsTd + n^UdW)^^ 



P(0„Od) 



\Z\\i - \nd\^ - \Td\^) 

+\nd\\\ - \Ti,\^ - \%\^) 

-(2|7^rf|2|7;|2 + T;T;n,nd + ZTdKiz*)!, 

(1 - \n,\^ ^ - iTZdl^ ~ \Td\^) 
HVUs + TV.TsXKTd + VTd)!, 



where we have introduced the (real) overlap integral 
dcj I dco'il/{cj,io')il/*{cj',io). 

-oo ^ — CO 
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